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Abstract 



One-dimensional three-body Schrodinger equation is studied, with binding mediated 
by the power-law two-body long-range attraction V^^\x) = a^x^ in superposition 
with a short-range repulsion V'^'^^ix) = P'^/x^ (plus, possibly, further subdominant 
power-law components). Such an unsolvable and non-separable generalization of 
Calogero model (where one had L = K = 2) is shown to become separable and 
solvable at all integers > and L > in the limit of the large repulsion (3^1. 



1 Introduction 



Solutions of an A— particle Sclirodinger bound-state problem 

^*(f 1, f2, ...,xa) = E'^{xi, X2,..., xa) (1) 

are vital in atomic and nuclear physics as well as in quantum chemistry. Using some 
realistic input Hamiltonian H one often calculates the low-lying spectrum of energies 
as a characteristic output. Its construction is not easy in general. Even the results 
of some very advanced (usually, variational) numerical methods need not be reliable 
and may require an independent verification, say, via an at least partially exactly 
solvable, simplified A— body model (their concise review appeared recently in [1]). 

Of course, the most natural opportunity for the (partial) exact solvability occurs 
in many-body Schrodinger equation (1) with the first nontrivial choice oi A = 3. In 
the single dimension with Xj e M and for the spinless and equal-mass particles, the 
solvability may play an independent important role in phenomenological as well as 
purely mathematical applications of the pertaining three-body models 

H^{xi,X2,X3) ^ E^{xi,X2,X3) H ^ -J^ + W{xi, X2, X3) , (2) 

i=l 

especially when their binding is mediated by the mere two-body forces, 

W{xi,X2, X3) = V{xi - X2) + V{X2 - X3) + V{X3 - Xi) . (3) 

In such a setting, in 1969 [2], Francesco Calogero discovered the exact solvability 
of the model which was based on the use of the most elementary spiked harmonic- 
oscillator (SHO) two-body forces 

W^''°)fe-x,)=a,^x<-x,)^+^^<^. (4) 

yxi Xj) 

This choice was motivated not only by the strong immediate phenomenological ap- 
peal of the SHO two-body potential itself but also by the encouraging fact that at 
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the vanishing z/ = 0, the simpUficd Schrodinger equation (2) - (4) is well known 
to remain separable. This observation opened the path towards the closed u ^ 
solution since in the polar Jacobi coordinates, its explicit construction emerges im- 
mediately after the reduction of the singular part of the sum (3) of potentials (4) to 
an elementary expression [2]. 

Our present paper may be perceived as a direct continuation of the same A = 3 
project, paying attention just to a broader class of potentials characterized by their 
more general power-law behaviour in both the short- and long-range domains. These 
potentials may be chosen in the similar two-term form 

V{xi-Xj)^a'^{xi-Xj)^ + P^{xi-Xj)'^, L,X = 1,2,... (5) 

or in the more general polynomial form (27) of Appendix A. They may mimic many 
popular confining forces, say, of the Lennard-Jones type. They combine a suit- 
able strongly repulsive singular core in the origin [3] with the standard anharmonic- 
oscillator asympotics. One arrives, in this way, at the most natural generalization of 
the Calogero's forces, indeed. 

We shall analyze the latter three-body problem in its two exphcit partial differ- 
ential forms (28) and (30) of Appendix A. The main purpose of such a study was, 
originally, a search for all the possible remnants of the exact solvability which is, 
definitely, lost in its completeness at any L > 2 and/or K > 2. The project found 
an initial encouragement also in our unexpected observation that the transition to 
L — A need not necessarily mean the loss of separability of eq. (2) (see section 2 
and, in particular, subsection 2.1 below). In parallel, we believed that even after 
a complete loss of separability, the problem might remain approximately separable 
and/or solvable. 

In this direction, two main technical ingredients of our study of eq. (2) consist in 
the use of the Jacobi coordinates at any K and L (cf . Appendix A) and in a successful 
generalization of the above-mentioned Calogero's trigonometric K — 2 identity to 
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all the integer exponents K and L (cf. Appendix B). An ultimate incentive for the 
presentation of all these results came from our a posteriori observation that the latter 
identities prove amazingly elementary and compact. 

In a way guided by a few examples revealing the simplifying role of the growth of 
P (and described in more detail in subsections 2.2 and 2.3) we are going to employ the 
simplicity of our auxiliary trigonometric identities of Appendix B for conversion or 
re-arrangement of our partial differential (PD) Schrodinger eq. (2) into a separable 
problem [i.e., into a set of two ordinary differential (OD) equations] plus corrections. 

This idea in its application to our non-separable ^4 = 3 problem is fully developed 
in section 3. We demonstrate there how the practical relevance of the separability- 
violating corrections decreases very quickly with an increase of the dominant repulsive 
component in superpositions (5). We show how our separable large—/? approximants 
may be constructed as solvable and how they generate the reliable low-lying three- 
body spectra. For an explicit illustration of the efficiency of such a strategy, we 
first test it on the solvable K = L = 2 case (i.e., on the model of Calogero) in 
subsection 3.1. 

In subsection 3.2 we move to the first "unsolvable" example with L — 3 and show 
how the comparatively compact character of the trigonometric identities of Appendix 
B enable us to achieve a full quantitative understanding of the competing mechanisms 
of repulsion and confinement in our non-separable two-dimensional effective potential 
wells with complicated shapes. 

Our final illustrative example of subsection 3.3 employs K — L — 3 and repre- 
sents, therefore, a full fledged modification of the Calogero's model in both its short- 
and long-range interaction parts. Of course, this modification remains non-separable 
and is also exactly solvable only in the limit (3 ^ oo. 

Section 4 is the summary of our results which re-emphasizes that all our three- 
body spiked anharmonic oscillators share the features demonstrated via particular 
examples. Thus, for all of them the low-lying spectra become available in the physi- 
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cally very appealing asymptotic domain characterized by the very strong short-range 
repulsion in our two-body interactions. 

2 Elementary guide: quart ic oscillator 

2.1 Quartic Schrodinger equation and its "forgotten" sepa- 
rability 

We need not recall the techiques presented in Appendices A and B to see that the 
potential function (3) remains virtually trivial, in the cartesian coordinates X and 
y, for all the three-body regular harmonic oscillators (RHO) with L = 2 and /3 = 
in (5). It is well known that the related PD Schrodinger equation (28) with 

u^'^\x,Y) = 3x^ + 3 

degenerates to a pair of the ordinary and solvable confluent hypergeometric differen- 
tial equations. In polar coordinates, in addition, this two-dimensional RHO potential 
term becomes completely independent of the angular variable </?, 

Still, to our great surprise we revealed that this example is not unique. There exists 
another two-body potential which leaves our PD Schrodinger equation (30) separable. 
Such an assertion can be easily verified since in the context of Appendix B.l we may 
derive that 

c/(^)(x,y) = ^(x^ + rf , 

i.e., 

n('\g,ip) = U('^[X{g,ip),Y{g,ip)] = ^g\ 

We see that the quartic oscillator belongs to a very exceptional set of three-body 
problems which are rigorously reducible to the solutions of an ordinary differential 
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equation. Such an observation will definitely deserve a separate treatment, say, in 
the context of perturbation theory. 

In what follows we are going to skip this opportunity because the separability of 
the quartic three body oscillator seems so exceptional and un-paralleled by any other 
power-law model. Indeed, in the light of the discussion in Appendix B we may be 
certain that the spiked quartic anharmonic oscillator (SQAO) two-body interaction 

V^'^^°\x, - xj) = u;' {x, - xjf + \{x,- xjf + f^^^ (7) 

[Xi Xj) 

generates the only set of the PD Schrodinger equations (30) which remain separable 
in the polar Jacobi coordinates. 

2.2 The role of the pronounced two-body repulsion 

In a marginal remark let us add that in contrast to the above-mentioned Calogerian 
K — L — 2 model, the separability of the latter quartic model does not imply its 
exact solvabihty of course. Although its two-body input force (7) is merely shghtly 
more complicated than its Calogerian SHO predecessor (4), the SQAO bound state 
problem must be solved in terms of perturbation expansions or by a numerical, 
approximate method [4]. 

In spite of the (expected) loss of the exact solvability, the spiked quartic illustra- 
tion may still profit from the above-mentioned experience gained during the study of 
the exceptional force (4). One learns, in particular, that the apparently complicated 
expression 

emerging in the Calogero's paper [2] may be significantly simplified by its further 
re-parametrization [cf. eq. (29) in Appendix A]. As long as we work now with the 
same repulsive core, the same function occurs in SQAO case. We have to return to 
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the Calogero's "experimentally" discovered trigonometric identity 

1 1 1 

+ —-1 ^7-7T + 



1/1 1 1 ^ _ 9 

2^2 l^sin^v^ sin^(v5 + |7r) sin^(v9 - |7r)/ 2^231,^23^ 



which makes, obviously, also the SQAO version of the PD equation (28) manifestly 
separable. 

Our fascination by the Calogero's model concerns, first of all, an ease with which 
one may replace the difficult, PD Schrodinger equation (2) by a system of independent 
OD equations. In this language, the introduction of the spherical Jacobi coordinates 
replaces both the SHO and SQAO versions of the PD problem (2) by the much more 
easily tractable "radial" OD equation 



^(r) = re (0,00) (9) 



with A = or A > 0, respectively. 

It is well known that we have here v — (^l — 3)/2 + £ where the integer index 
£ = 0, 1, . . . numbers the so called (hyper) spherical harmonics [5]. One also recollects, 
vaguely, that an addition of a new force Vspike^r) ^ gjr^ would remain tractable in 
eq. (9) because its form coincides precisely with the kinetic (so called centrifugal) 
singular term so that we just have to admit any real v — v{g) [3]. 

For the time being, let us fix A = and pay attention just to the exactly solvable 
spiked harmonic oscillators with a strong repulsion g ^ \ (i.e., with a large real 
V ^ \). Of course, their spectrum En = uj (An + 2^ + 3) with n = 0, 1, . . . is 
equidistant and moves merely upwards with the steady growth of i^. Even without 
any knowledge of the exact formulae, this phenomenon is easily explained by an 
elementary observation that near the minimum, the shape of the SHO potential term 
in eq. (9) may be very well approximated by a one-dimensional RHO well. Thus, we 
may starts from the formal Taylor series 

Veffir) = ^^^^ + a;V = V,ff{R) + ^(r - Rf V^f^iR) + 0[{r- Rf V^j^iR) 
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2u;^R^ + Au;'{r-Ry + 



(r - Rf 
R 



(10) 



In the vicinity of the absolute minimum at r = R = + l)/a;^]^/*^ and in the 
asymptotic domain characterized by the "thick" spike term, v 1, the latter po- 
tential is very well approximated by its first two terms since ~ l/v^ ^ 1- 

One arrives at a consistent picture where the approximate energy levels calculated 
for the one-dimensional leading-order RHO well, i.e., the values 

E^^Veff{R)+2uj{2m + l), m = 0, 1, . . . 

coincide with the above exact low- lying SHO spectrum En,Ei,... up to the second 
order term in 1/R. This is an extremely encouraging observation which shows that 
the domain where the short-range repulsion is large is, in a way, privileged. 

One should keep in mind that in our exactly solvable illustration, the parameter 
1/R is precisely the appropriate measure of smallness of corrections. Once its value 
proves sufficiently small, we arrive at an asymptotic equivalence of the two models. 



RHO 



SHO. 



(11) 



In addition, we may easily move beyond the leading-order precision, making use of 
the textbook perturbation series [6]. Their explicit construction to arbitrary order 
and for virtually any potential may proceed, e.g., along the lines described in full 
detail in the review paper [7]. 

Same observations may be easily extended to the "purely numerical" SQAO prob- 
lem where the effective potential Veff{r) — ^^^^r^ -|- a;^r^ + Xr^ with r — R + x may 
equally well be expanded in the similar Taylor series 



+ cu^R^ + XR' 



i?2 



+ 2 



u^R ^—^ + 2XR^ 



x + ... . 



We see that the point R becomes an absolute minimum of Veff{r) whenever the 
linear, 0{x) term vanishes. This specifies the value of R exactly (i.e., as a positive 
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root of the cubic algebraic equation) and enables us to find the subsequent quadratic 
and higher Taylor-series terms 



+ 



+ XR 



X 



+ 0{x' 



12) 



in closed form. 



2.3 Another example possessing an inverse-cubic core 

In previous section wc have chosen the exactly solvable illustration just to emphasize 
the key idea of construction which lies in the coincidence of the shapes of two different 
potentials near their minima, caused by the growing strength of the singularity in 
the origin. Let us now verify an applicability of such an approach using a slightly 
less trivial form of the single-particle repulsion, 

^(r) = 0, W{r)^Fr^ + ^, r, F, G > . (13) 

At a large value of G » 1, potential W{r) may be represented by its Taylor series 
again. Near its absolute minimum which occurs at rmin = R = (G/F)^/^ » 1 we 
have 

W{R + ^)^2^/GF + 9^/GF^-9^/GF-^ + ... . (14) 

This expansion illustrates the feasibility of the same harmonic-oscillator approxima- 
tion approach to our generic "unsolvable example" (UE). In the light of the large— z/ 
equivalence SHO RHO between the two exceptional but "equally exactly" solv- 
able examples, the unsolvable equation (13) and its asymptotic solvability opens a 
way towards a triple asymptotic equivalence between the three low lying spectra, 

UE < — > RHO < — ^ SHO. (15) 

This means that we may, alternatively, approximate also any unsolvable model by 
the zero-order solvable model with a singularity. Of course, the triple-osculation 
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property (15) restricts our choice of the spiked SHO model characterized, for clarity, 
by the zero subscript o, 

W^o(i?o + = ^o(i?o + 0' + 7^-^^-2yG^ + 4Foe'-4(/^e' + --- ■ (16) 

[Ro + V <^o 

We must require that the minima are the same, Rq = R, and that the two harmonic 
oscillator parts of the wells (14) and (16) differ at most by a convenient constant 
shift. These requirements form a set of two equations with easy solution. 

In the domain of the large G ^ 1 we have Go ^ -^o as expected. 

We may summarize that all the unsolvable OD examples may be studied via 
different implementations of the same idea of the strong-repulsion approximation. 
Of course, even when all the available free parameters coincide as they should, one 
may still test and compare the practical numerical performance of the RHO and 
SHO alternatives. This will not be pursued here. We only shortly emphasize that 
the use of the more sophisticated two-parametric SHO zero-order forces with a core 
gives a better chance for an optimal choice of the zero-order approximant. A priori, 
an additional benefit might be spotted in the better coincidence of the domains 
(i.e., the half-axes of the coordinates) which might also support the preference of 
the pair of the UE and SHO Hamiltonian operators. Thus, the similarity of the 
shape (often called osculation) of the given curve W'^^^^ near its minimum with the 
solvable spiked well Wq^^^^ appears to be better founded, in spite of its shghtly 
less immediate construction. Presumably, the importance of the similar arguments 
might further increase at ^4 = 3, in the genuine PD context. 
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3 Three-body spectra at any K and L 



3.1 Solvable guide: The K = L = 2 model at z/ > 1 

The only non-constant angular dependence in the Calogero's Q^^^'^\g, (p) occurs in 
its short-range repulsion component Q^~'^\g,(p) ~ sin ~^ 89?. For ip E (0,7r/3) (i.e., 
within a physical wedge in the phase space which is called, sometimes, the Weyl's 
chamber), this function has a unique minimum at (fimin = 7r/6 where sin3</7^i„ = 1. 
In its vicinity, the angular Taylor series may be calculated, 

Sm 3(7 -|- TT/Dj 5 

This means that the absolute minimum of the potential fl^^"°\g,(fi) will he at 
^ — Vmin = 7r/6. For the determination of its second coordinate g — R, we may 
analyze the function Q'^^^'^\g,TT/Q) = of the single variable ^ > 0. The 

second necessary Taylor series is then very easily derived, 

Veffie) = ^^^^^ + ^^'q' = VeffiR) + \{q' Rf v:,,{K) + o[{g- Kf v:if{R) 

g-Rf ' 



6u;'R' + 12u;'{g- RY + O 



R 

Obviously, we just have to assume that i? 1 in order to achieve a full consistency 
of the Calogerian model with its two-dimensional RHO separable approximation, 

n^^^^\R + ^,7]/R + 7r/6) ^ Guj'^R^ + 12(jj^^'^ + 27 uj^ 7]"^ . (18) 

In a way predicted by the single-particle guide of section 2 one discovers that all the 
higher-order terms 0{ri"^ x ^"') remain suppressed by the same factor <^ i 

Thus, in the domain of large u, the low-lying spectrum of energies may be approxi- 
mated by the RHO formula 

E = Em,n = Q^^^R^ + ^VSou {2n + 1) + 3V3u}{2m + 1) + . . . . (19) 
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Up to the higher-order corrections, this prediction "explains" the well known exact 
spectrum of the low- lying energies in the Calogero model [2] . 

We may conclude that in spite of the use of the polar coordinates, the model 
is suitable for the two-dimensional RHO approximation using the original cartesian 
coordinates rotated merely by 7r/6. This is an immediate consequence of the structure 
of the kinetic term in eq. (30), 

_ 92 _ 1 

Obviously, it may be locally re-interpreted as a cartesian kinetic energy, 

with the re-scaling of one of the coordinates which we employed also in eq. (18) 
above. 

3.2 A = S confinement and its sample proof for the non- 
separable L = 3 oscillator 

One of the simplest non-separable three-body models may be based on the cubic 
anharmonic-oscillator two-body potentials 

V(^^^^\x, - xj) = {x, - xjf + 7 (^. - xjf + f^^^ . (20) 

[Xi Xj) 

Of course, one must proceed with a certain care since even in the single-particle 
case, the cubic potential is not confining on the real line [8]. In the present setting, 
fortunately, we need not work on the whole real line. On the contrary, due to the 
presence of the strong repulsion, we must parallel the Calogero's considerations [2] 
and fix the ordering or numbering of our three particles in advance, keeping in mind 
that in quantum case, our particles cannot tunnel through the l/{xi — XjY barriers. 
This means that we may demand that, say, xi > X2 > X3. This automatically 
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guarantees that our system will "live" in a single "physical" wedge with G (0, 7r/3). 
At this moment we may also guarantee that the cubic component fl^^\g,ip) of the 
whole potential Q^^^^'^\g, </?) will remain non- negative (i.e., confining or vanishing) 
at all the positive couplings 7 > 0. 

In a certain challenging extreme, let us now drop the safely confining quadratic 
force in eq. (20) and show that even at c<j = 0, the thrcc-body bound states will 
still remain localized. In essence, we have to demonstrate that the spiked cubic (SC) 
two-body forces V^^'^\xi — Xj) [— V'^'^'^^'^^Xi — Xj) at a; = 0] will guarantee the 
confinement of the three-body system even at the ends of the segment (p e (0, 7r/3) 
where the potential Q^^\g,(p) itself vanishes. 

The proof proceeds as follows. Firstly, assuming that g is sufficiently large we 
show that with respect to the angle if, there exist two minima in the complete SC 
potential 

n^''^\g,<p)^a'g'sm3ip+ ^ ^\ , > ^ ^^^+1) > . (21) 

^^sm 6ip 1 I 

They occur at the two ^—dependent angles </7± = 7r/6 — </7=p. Their ^—dependence is 
controlled by the rule d^p^'^^'^^ — with the following solution and its consequence, 

52 \ 1/5 



sin^3^± = ^ ^ 9.^^^)>^{2a'i5-g)"\ g>g,= [^ 



2^ 
a2 



This means that the potential ^^^'^^ grows sufficiently quickly as a function of g in 
the whole g^\ asymptotic wedge. The model possesses just the discrete spectrum. 

Marginally we may note that the latter proof may be generalized from the above 
very special cubic model to all the potentials r2*^^*^+^^(f), (/?). Indeed, as long as the 
angular minimum moves quite quickly to the wedge boundary with the growing M , 
</?_ = 7r/3 — ~ ^-(2A^+3)/3^ ^jjg insertion of this value in any M > 1 potential 
produces the general estimate 
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We see that for all the integer exponents M > 0, no particle can escape in infinity. 
The potential grows at least as a 4M/3— th power of q in the whole interval of the 
angles in the asymptotic domain of q. 



3.3 An ultimate illustration: A = 2> spectrum at = L = 3 

The climax of our analysis and considerations comes with the strong-repulsion con- 
struction of the energies in any power-law A — ?> system. For the sake of brevity 

we shall pick up, pars pro fMo, the same unsolvablc spiked cubic oscillator as above, 
described by the non-separable PD Schrodinger equation 



We already know that its potential has a double-well shape with two minima at 

every fixed and sufficiently large radial distance g ^ 1 from the origin. With the 
decrease of g, both these minima decrease and their positions move towards the 
central maximum (which "sits" at (p — (po — n/G). All these three extremes merge 
at ^ = ^0 = Below this distance, just a single minimum survives in the 

middle of the valley, at cp = (po. 

The potential at its minimum still decreases with the further decrease of g, and 
the unique absolute minimum of ^l^^'~^\g, cp) is finally reached at (p = (po = 7r/6 and 
at the distance 



This quantity is large whenever ^ 1 so that we may use it as a new, slightly 
simpler measure of the strength of the SC two-body repulsion. This means that we 
may re-define = 3 g — R + ^, (p — 7r/6 + rj/R and expand 



Q2 1 Q2 



- + n'^sc){g,p) ij{g,^) = Eij{g,ip), (22) 



dg'^ dip' 



R = ^("^^") = ^2/3V(3q;2). 




o? g^ siv? 3(/? 
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We see that the corrections decrease with the integer powers of 1/ R. This enables 
us to insert their scries in eq. (22). With a new re-scahng constant a = l/\fR 
introduced in the two independent variables, 



Q = R + a X, 



= 7r/6 + ay/R, 



this gives our final approximate Schrodinger equation 



dx^ dy^ 



(23) 



The reliability of this approximation will increase with the growth oi R ^ 1 or 
a = a{R) < 1. 

At this point one must emphasize the two parallel merits of the approximation 
(23). Firstly, its is a separable PD equation which may be analyzed as two indepen- 
dent OD equations. Secondly, both these OD components are exactly solvable. This 
means that in the strongly spiked domain, the low-lying energies are given by the 
following closed formula 



We see that these energies are numbered by the respective "radial" and "angular" 
quantum numbers m, n = 0, 1, . . . in a way which resembles the previous Calogerian 
exactly solvable case. Thus, one may expect that the same asymptotic construction 
remains applicable to all our models (28). 



m,n 




4 Summary 



By our present paper, two main messages are delivered: 
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• the power-law confinement of tfiree spinless particles in one dimension leads to 
a comparatively simple PD Schrodinger equation, provided only that we work 
in the spherical Jacobi coordinates; 

• a quantum osculation method generates the low- lying spectra with the precision 
which increases with the strength of the repulsive two-body core in V{xi — Xj). 

An exceptional, unexpected separability of the quartic anharmonic oscillations of 
three particles has been also revealed as a byproduct of our systematic explicit con- 
struction of the corresponding PD Hamiltonian operators. 

The former pair of messages resulted, basically, from a successful generalization 
of the Calogero's trigonometric identity (33). In the historical perspective of ref. [2], 
the latter identity opened the path towards the explicit construction of the whole 
class of solvable models. Some of the consequences of its present generahzations are 
equally exciting: This is the essence of our mathematical achievement. 

One of its most important uses for physics may be seen in the context of the 
large— /3 expansions. Their construction at any K ^ 2 ^ L seems to be significantly 
facilitated by a generic survival of the separability beyond the Calogero's exceptional 
model whenever we move to the strongly repulsive regime. In this sense, even all the 
multi-term generalizations of our present power-law A = ?> oscillators characterized, 
in both the long-range and short-range regime, by a polynomial behaviour of their 
two-body potentials remain tractable as models which will be approximately separable 
and solvable in the domain of the very strong two-body repulsion at short distances. 

One of the most remarkable purely technical aspects of our present work appeared 
to be an unexpectedly slow growth of the complexity of the closed trigonometric 
formulae for the total potentials n^"^\g,(f ) with the increase of absolute values of 
their maximal-power superscripts m. This is a lucky circumstance which enhances 
the feasibility of the construction of the approximate spectra in polar coordinates 
and in the strongly-spiked limit significantly. 
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In the future, one may expect that also a systematic incorporation of the higher- 
order corrections will remain feasible. We plan the more systematic study of such 
a possibility based on a refined strong-coupling perturbation-series expansion of the 
observable quantities in the powers of the auxiliary quantity 1/R where R denotes 
the radial distance of the absolute minimum of fi("^)(^,<^) and grows in proportion 
to the strength of the repulsive two-body core. 
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Appendix A: Jacobi coordinates at A = 3 



It is well known that the introduction of the suitably normalized centre-of-mass 
coordinate Z = Z^"^-* = Yl,k=\ ^kl^~^ enables the elimination of the free bulk motion. 
Of course, the use of the new variable Z imposes a constraint upon the old >1— plet 
of coordinates. The redundancy of one of the positions may be easily resolved by 
the transition to the multiplet (Z, X, Y, . . .) of the so called Jacobi coordinates which 
may be defined at any number of particles A and which we shall specify, at A = 3, 
by their most transparent matrix definition 

^ 1/V3 ^ 

l/\/2 -l/\/2 

y 1/a/6 1/V6 -2/a/6 ) 

with an elementary inversion. In this manner, the original particle coordinates may 
be understood as a triplet of linear functions of three new parameters, 

/ 



X 



1.^ 



(25) 



xi{Z,X,Y) 
X2{Z,X,Y) 
Xs{Z,X,Y) J 



^ ( 1/V3 1/V2 1/VG \ ( Z^ 



\ 



J 



X 



(26) 



1/V3 -1/^2 l/\/6 

1/a/3 -2/^/6 
These functions enter the general interaction term in eq. (2), 

L 

W{xi,X2,xs)= J2 FmW^"'\xi,X2,xs), = > 0, F_K^P^>0, 

m=-K 

W^'^^X,, X2, X3) = {Xi - X2r + (^2 - Xsr + (^3 - X^r (27) 

where the positivity of guarantees the confinement while the positivity of F_k 
prevents the system from a collapse. In such a notation we may re-write our ^4 = 3 
Schrodinger equation as a reduced PD problem in the two cartesian coordinates X 
and y. 



+ J2 FmU^"'\X,Y)-E \ ^X,Y) = 



(28) 



m=-K 
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where the functions 

U^^^ {X, Y) = W^^^ [x,{Z, X, Y),X2{Z, X, Y),x,{Z, X, Y)] 

prove independent of Z. One may also employ the polar re-parametrization of the 
Jacobi coordinates, 

X — X{g,(p) — g sm(p, Y — Y{g,ip) — g cos cp , (29) 

arriving at an alternative, equivalent formulation of the same Schrodinger equation, 

dg^ ^2 Q^2 



+ F^n^'^^\g,^)-E\ ^X{g,^),Y{g,^)] = (30) 



where we introduced another abbreviation 

n^^\g,^) = U^^^[X(g,^),Y(g,^)]. 

Appendix B. Auxiliary trigonometric identities 

Any practical analysis of the Schrodinger A — 3 equation (2) in Jacobi coordinates 
requires the explicit specification of its interaction term (3). Separately, let us derive 
its form for the three specific subclasses of the underlying two-body input force. 

B.l. Potentials t/^^^ and Q^"') with even m = 2M 

In our present notation, the regular two-body interactions are characterized by the 
even integer superscripts and by the following asymptotically confining interaction 
functions, 

C/(2^)(X,y) = [x^{R,X,Y)-X2{R,X,Yf + 

+ [X2{R, X, Y) - x^iR, X, Y)f^' + [xs{R, X, Y) - x^{R, X, Y)f^ . 

The growth of M means a strengthening of the confinement in both infinities, 
which has a definite phenomenological appeal. The explicit form of the functions 
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U^'^^\X,Y) is, unfortunately, more and more complicated for the larger and larger 
M, 

Y)^-X' + - X^Y' + — X'Y^ + - y^ 
4 4 4 4 

[/(^) {X, y) = ^ (^^ + (27 Y^ + 225 X^Y^ - 15 X + 43 X^) , 

etc. Still, our search for the parallels between the solvable and unsolvable models 
finds its reward here since with the polar coordinates where we introduce an abbre- 
viation 

= {xi [R, X{g, V9), Yig, ip)] - [R, X{q, 99), Y{q, + . . . 

the fairly perceivable simplifications emerge after the patient trigonometric manipu- 
lations. This is one of our most important technical results. We get the sequence of 

the pleasantly compact formulae, 

t 

4 V- ■ 4 



Q(6) = ^ (9 + 2 sin^ 3(^) = ^(10 - cos 6 (^), 
99) = ^ (27 + 16 sin^ 3<^) , 



97 nlO 

Q(io)(^^<^) = fi^ (9 + 10 sin2 3(^) 



16 

o(i2)r ^ 729^ 81^12 _ 3 12 _ 

3^ 3(^+^pSm 3(^, 

etc. A remarkable regularity emerges in this pattern. For example, the choice of the 
superscripts m — 2M — 6N, 6N -|- 2 or 6N + 4 leads to the compact general formula 

= ^ + sin^3(^ + ... + c^sin2^3(^ (31) 

with a independence on the right hand side and with an extremely elementary 
N— and M— dependence of some its coefficients which may be derived very easily 
whenever needed. 
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B.2. Potentials U^""^ and Q^"^) with odd m = 2M + 1 

Systematic study of all the potentials U^"^^ with the odd powers m = 3, 5, ... is not 
so easily extended to A = 4 but in the present A = 3 setting it is well defined and 
worth the study. The elementary experimenting teaches us quickly that up to the 
trivial U^^\X,Y) = we have to deal with the interesting functions. One cannot 
be surprised by the observation that the complexity of the representation of the 
complete potentials in Jacobi coordinates grows quite quickly with their superscript, 

c/(^) (X, y ) = ^ V2 (x^ - 3 y 2) + X, 

U^'\X, Y)^^V2(X^- 3y2) X (x^ + Y^)\ 

{X, Y)^^V2 (81 Y^ + 279 X^Y^ + 219 X^Y^ + 85 X^) (x^ - 3 F^) X 

etc. Fortunately, in the same manner as above, the use of the polar coordinates 
makes the geometric structure and symmetries of these two-dimensional potential 
wells much more evident, 

2 

n5 



n^^\g,cp) = 1^ V2 sin3(^, 
8 



3 0^ 

q(9) ip) = V2 sin 3(p (81 + 4 sin^ 3(p) 
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3 

V2 sin 3ip (83 — 2 cos 6 ip) 



33 o^^ 

^(11) ip) = V2 sin 3^ (27 + 4 sin^ 3^) , 



117 o^^ 

(^) = 2 ^ sin 393 (27 + 8 sin^ 3<^) , 



64 
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etc. Again, once we have m = 6N + 3 or m = 6A^ + 5 or m = 6A^ + 7, the exphcit 
form of the potential O^"^) will be prescribed by the trigonometric formula 



2M+r = + 1) \/2 sin i^^ + di sm^ 3^ + ... + dN sin^^ 3<^ 



where, formally. 



2M 



(32) 



N — entier 



M -1 



M = 1,2, 



and entier [x] denotes the integer part of a real number x. 



B.3. Strongly singular potentials Q^^^ with negative 



m 



Using the same abbreviations as above we may immediately complement the 
Calogero's identity (8) by its Coulombic predecessor, 

3 



( 3 



Q \f2 sin 393 ' 
\ 2 



(33) 



\ Q -\/2 sin 3<^ / 

These two formulae exhaust the set of the standard singular forces. Nevertheless, 
there is no physical reason for avoiding the more singular repulsion near the origin, 
and the simplicity of the latter two formulae attracts attention to the more general 
spikes in 

n(-'=)(^, ^) - (xi - X2)-' + (X2 - X3)-' + (X3 - xi)-'^ . (34) 

At the first three "nontrivial" exponents A; > 2 we reveal again the same overall 
structure of the interaction terms in their trigonometric representation, 

Q(-^)(^,^) = \p.^-'\QM'\ (\-\ Sin' 3^) , 
n(-'\g,ip) = [n^-'\g,ip)]' (1 - ^ sin^ 3 9.) , 
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n(-'\g,^)^[n(-'\e,^)f (i-^^ sin's 

Obviously, the pattern remains precisely the same as above. In place of any routine 
classifications, let us only add the first less trivial formula 



904 



824 



sin^ 3 99 + sin'^ 3 99 ) , 



667 18009 

which nicely illustrates not only the emergence of the new terms but also the close 
parallel between the positive and negative superscripts. 
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